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SECTION  I 


INTRODUCTION 


The  design  of  equipment  to  withstand  dynamic  loading  of  its 
parent  structure  has  been  addressed  for  a wide  range  of  structures  by 
means  of  theory  and  experiment.  Shock  effects  on  equipment  in  hardened, 
protective  structures  have  been  estimated  by  constructing  envelopes  to  the 
spectral  response  at  the  attachment  points,  Reference  1.  It  is  usual  to 
avoid  mounting  sensitive  equipment  directly  on  the  structure  if  the  fixed 
base  frequency  of  the  equipment  falls  near  a natural  frequency  of  the  structure; 
potentially  damaging  overlap  is  avoided  by  redesigning  the  equipment  or  by 
using  shock  isolation  systems.  References  1 and  2.  In  order  to  study  the 
effectiveness  of  shock  isolation  systems  in  complex  structure  over  a wide 
frequency  range,  impedance  testing  procedures  have  been  developed,  References 
3 and  4. 


I 
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Problems  involving  equipment-structure  interaction  have  been 
addressed  in  submarines,  where  the  equipment  is  massive,  and  may  be 
heavier  and  have  lower  natural  frequencies  than  the  structure,  References  5 
and  6.  Special  analysis  and  design  procedures  have  been  developed  to  take 
advantage  of  the  unique  characteristics  of  such  systems.  The  spectral 
response  at  the  attachment  point  of  massive  equipment  exhibits  a minimum  at 
its  fixed  base  frequencies.  Therefore  the  envelope  to  the  spectral  response 
is  unnecessarily  conservative  and  the  maximum  stresses  in  the  equipment 
will  most  probably  occur  at  the  minima  in  the  shock  spectra;  this  is  taken 
into  account  in  design,  Reference  7. 

Earthquake  response  of  structures  and  their  contents  is  another 
area  where  equipment-structure  interaction  has  been  addressed.  One  common 
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approach  is  the  floor  spectrum  method,  in  which  interaction  is  completely 
neglected.  The  equipment  is  treated  as  a single-degree-of-f reedom  system 
subject  to  the  base  motion  that  the  structure  would  experience  at  the 
attachment  point  in  the  absence  of  the  equipment.  In  addition  to  neglecting 
interaction,  time  history  analysis  of  the  structure  to  determine  that  motion 
must  be  carried  out.  Approximate  techniques  have  been  proposed  to  develop 
floor  response  spectra  from  ground  spectra  that  bypass  the  computational 
problems  associated  with  this,  but  these  are  ad  hoc  methods  and  no  estimate 
of  their  accuracy  can  be  made  (Reference  8) . An  alternative  approach  is  to 
consider  a N + 1 degree-of-freedom  model  for  the  equipment-structure  system 
and  subject  it  to  time  history  analyses  for  a variety  of  specified  ground 
motion  inputs.  This  can  be  cumbersome  and  expensive  for  design.  Further- 
more, conventional  methods  of  dynamic  analysis  are  accurate  in  the  computation 
of  the  response  for  the  lower  modes,  and  they  frequently  possess  some  form 
of  numerical  dissipation  to  damp  out  spurious  participation  from  the  higher 
modes.  In  the  cases  where  structure-equipment  interaction  is  important,  at 
intermediate  and  high  frequencies,  the  use  of  such  codes  can  mask  significant 
equipment  response.  Finally,  response  spectra  methods  for  N + 1 degree-of- 
freedom  systems  suffer  from  uncertainty  in  how  to  combine  the  peak  modal 
values.  The  approach  using  the  square  root  of  the  sum  of  the  squares. 
References  9 and  10,  has  been  shown  to  be  inaccurate  in  the  case  of  closely- 
spaced  modal  frequencies  and  light  equipment  having  a natural  frequency  equal 
to  one  of  the  natural  frequencies  of  the  structure.  Penzien  and  Chopra 
(Reference  11)  and  Penzien  (Reference  12)  have  studied  this  problem  for  a 
simple  tuned  system  and  have  proposed  a method  in  which  the  response  spectra 
of  two-degree-of-f reedom  systems  are  used.  In  this  approach,  the  N modes  of 
of  the  structure  are  determined  and  to  each  mode  considered  as  a single- 
degree-of-freedom  system,  the  equipment  is  attached  to  give  a two-degree- 
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of-freedom  system.  The  time  history  of  the  response  of  each  of  these  N 
two  degree  of  freedom  systems  is  computed  numerically  and  the  maximum 
response  of  the  entire  system  evaluated  using  the  square  root  of  the  sum 
of  the  squares  of  the  maximum  value  for  each.  This  method  is  shown  by 
comparison  with  a direct  computation  of  the  time  history  of  the  response 
of  the  total  system  to  be  fairly  accurate.  It  is  shown  below  that  the 
method  is  unnecessarily  complicated,  in  principle  incorrect,  and  that  its 
accuracy,  established  in  Reference  11  only  for  the  particular  low  order 
system  studied,  is  in  a sense  fortuitous. 

The  present  work  is  concerned  with  a rational  approach  to  one 
aspect  of  this  problem,  that  of  light  equipment  whose  frequency  is  close 
to  one  of  the  natural  frequencies  of  the  structure,  a situation  described 
as  "tuning."  To  gain  insight  into  this  problem  area,  idealized  models 
are  considered  which  incorporate  the  characteristics  of  a structure  and 
internal  equipment,  the  equipment  having  natural  frequencies  which  are 
higher  than  the  fundamental  frequency  of  the  structure.  The  combined 
structure-equipment  system  considered  is  one  in  which  the  equipment  is 
relatively  light.  In  previous  work  (Reference  13),  we  have  described  the 
response  of  this  system  to  steady-state  ground  shaking.  It  was  shown  that 
significant  interaction  effects  occur  only  in  the  case  characterized  as 
tuning,  or  near-tuning;  namely,  the  situation  where  the  frequency  of  the 
equipment  considered  as  a single-degree-of-f reedom  system  is  the  same  as, 
or  close  to,  one  of  the  natural  frequencies  of  the  structure.  If  the 
equipment  frequency  is  not  tuned,  or  nearly  tuned,  to  a structural 
frequency,  the  response  curve  is  roughly  the  superposition  of  the  structural 
response  and  the  equipment  response  with  little  interaction.  This  indicates 
that  the  conventional  floor  spectrum  method  is  valid  for  the  solution  of 
transient  problems. 
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If,  on  the  other  hand,  the  equipment  frequency  is  tuned  to  a struc- 
tural frequency,  it  was  found  that  for  the  combined  system  there  were  two 
closely  spaced  frequencies  on  either  side  of  the  tuning  frequency  and  that 
a band  of  high  amplification  appears  around  the  tuning  frequency  which  offers 
a substantial  target  for  sympathetic  oscillation.  Thus  there  is  a signi- 
ficant interaction  in  this  situation  between  the  structure  and  the  equipment 
and  suggests  that  the  conventional  floor  spectrum  method,  which  ignores 
such  interaction,  will  not  be  valid  for  the  transient  analysis  of  tuned,  or 
nearly  tuned,  systems. 

A typical  result  for  the  steady  state  response  of  a structure- 
equipment  system  is  shown  in  Figure  1 which  is  taken  from  Reference  13.  In 
this  case,  the  equipment  was  tuned  to  the  third  natural  frequency  of  the 
structure,  and  the  curve  shown  is  the  ratio  of  the  equipment  acceleration 
to  the  input  ground  acceleration  considered  as  a function  of  the  frequency 
of  the  input  normalized  with  respect  to  the  natural  frequency  of  the 
equipment.  The  mass  of  the  structure  in  this  example  was  a thousand  times 
that  of  the  equipment. 

In  this  report,  we  will  describe  a method  valid  for  the  transient 
analysis  of  tuned  or  nearly  tuned  structure-equipment  interaction.  This 
method  is  considered  as  a replacement  for  the  floor  spectrum  analysis.  It 
utilizes  the  shock  spectra  for  the  specified  input  to  the  structure  and 
directly  provides  the  shock  spectra  for  the  equipment.  In  this  sense,  it  is 
simpler  and  easier  to  apply  to  design  than  the  floor  spectrum  method  which 
requires  the  time  history  of  the  structure  response  to  be  computed  in  order 
to  obtain  the  equipment  spectrum.  The  simplicity  of  these  results  is  due 
to  the  fact  that  it  is  possible  to  take  advantage  of  the  mathematical 
structure  of  the  equations  of  the  tuned  or  nearly  tuned  system  and  also  to 
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make  use  of  asymptotic  methods  made  possible  by  the  smallness  of  the  equipment 
mass  in  comparison  to  the  structure  mass. 

The  system  to  be  studied  is  one  for  which  the  mass  of  the  equipment 
is  small  in  comparison  to  thr-  of  the  structure.  Damping  in  both  the  equip- 
ment and  the  structure  is  also  small.  The  equipment  is  taken  as  tuned,  or 
nearly  tuned,  to  a natural  frequency  of  the  structure.  It  is  assumed  through- 
out that  the  information  available  to  the  designer  is  a design  spectrum  for 
the  ground  motion  and  fixed-base  modal  properties  of  the  structure  alone  and 
of  the  equipment  alone.  The  results  obtained  will  be  an  estimate  of 
the  maximum  displacement  and  maximum  force  experienced  by  the  equipment. 

The  rationale  for  using  design  spectra  methods  is  that  they  are  cheap  and 
to  a certain  extent  include  the  probabilistic  nature  of  the  structure 
definition,  i.e.,  the  lack  of  precision  in  specifying  the  structural  parameters 
and  the  probabilistic  nature  of  shock  or  earthquake  input.  The  probabilistic 
nature  of  the  problem  is  taken  into  account  in  the  construction  of  a design 
spectrum  itself  and  in  the  way  the  maximum  values  in  each  mode  are  combined  to 
predict  the  maximum  for  the  entire  system. 

Under  the  limitations  outlined  above,  the  results  obtained  are  easy 
and  efficient  for  practical  use  by  a designer.  A feature  of  the  present 
analysis,  is  the  extremely  simple  result  that  if  response  spectrum  for  the 
ground  motion  is  available,  the  response  spectrum  for  the  equipment  can  be 
calculated  merely  by  multiplying  it  by  an  amplification  factor  which  depends 
on  the  mass  ratio,  the  degree  of  detuning  and  the  damping. 

The  reason  for  the  simplicity  of  the  result  can  be  indicated  on 
physical  grounds,  if  it  is  recognized  that  the  major  portion  of  the  response 
is  a classical  beat  phenomenon.  It  is  well  known  that  in  weakly  coupled 
systems  with  the  same  frequency,  the  response  of  the  system  involves  a per- 


here.  The  weak  coupling  is  achieved  through  the  small  mass  ratio  of  equip- 
ment to  structure.  For  simplicity  of  explanation,  consider  an  undamped  tuned 
equipment-structure  system.  Were  the  structure  alone  subjected  to  a ground 
motion,  the  velocity  imparted  to  it  would  be  independent  of  the  mass  and 
determined  only  by  the  ground  motion.  If  the  same  ground  motion  were  to  be 
applied  directly  to  the  equipment  alone,  the  same  velocity  would  be  imparted 
to  it  in  the  case  of  tuning.  The  kinetic  energy,  on  the  other  hand,  will  be 
proportional  to  the  mass  of  the  system  excited.  In  the  case  of  the  equipment, 
this  would  be  much  smaller  than  that  of  the  structure.  However,  if  the  equip- 
ment is  attached  to  the  structure  and  the  structure  is  subjected  to  ground 
motion,  the  kinetic  energy  imparted  to  the  structure  is  transmitted  in  its 
entirety  to  the  equipment  if  tuned,  and  the  velocity  imparted  is  amplified 
by  the  reciprocal  of  the  square  root  of  the  mass  ratio.  This  simple  result 
recurs  for  all  other  response  spectra  of  the  equipment. 

It  is  clear  that  damping  will  play  an  important  role  in  this  pro- 
cess, since  the  energy  transfer  takes  many  cycles  and  much  of  the  kinetic 
energy  in  a damped  system  could  be  dissipated  before  being  transferred. 

This  is  accounted  for  in  the  analysis  when  damping  is  included.  We  have  also 
studied  the  case  when  the  equipment  frequency  is  not  perfectly  tuned  to  a 
natural  frequency  of  the  structure.  If  the  detuning  is  small,  then  two 
closely  spaced  frequencies  again  occur  and  the  determination  of  the  resulting 
peak  response  must  be  calculated  with  care,  as  for  the  tuned  case.  The  final 
result  is  analogous  to  the  tuned  situation,  but  the  amplification  factor  is 
modified  by  the  inclusion  of  a parameter  describing  the  degree  of  detuning. 
When  the  equipment  frequency  is  grossly  detuned  from  that  of  all  natural 
frequencies  of  the  structure,  then  the  floor  spectrum  method  is  valid. 
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However,  we  have  derived  a much  simpler  method  which  gives  the  peak  response 
of  the  equipment  in  terms  of  the  given  ground  spectra  multiplied  by  amplifi- 
cation factors  depending  upon  the  fixed-base  modal  properties  of  the  structure 
alone  and  of  the  equipment  alone. 


SECTION  2 


SUMMARY  OF  RESULTS 


To  assist  the  reader  in  understanding  the  theoretical  development 

and  its  applications,  the  final  result  is  presented  in  this  section  in  terms 

of  the  free-field  spectral  acceleration,  S^,  which  is  a function  of  frequency 

and  damping,  and  a modification  factor,  which  depends  on  properties  of  the 

equipment  and  the  structure.  Although  the  theory  is  developed  for  the 

special  cases  of  systems  which  are  (a)  damped,  tuned;  (b)  undamped,  slightly 

detuned;  and  (c)  damped,  slightly  detuned,  the  result  is  given  by  one  general 

expression  as  follows: 

SA  ( (1+§)„. 

^ U 

(Y+S  +ABB)'1 


where 

* free-field  spectral  acceleration  evaluated  at  the  frequency 
(1  + ^)o>,  with  spectral  damping  ratio  (8+B)/2. 

£ = detuning  parameter,  expressing  the  degree  of  detuning 

between  the  structure  (mode  of  interest  has  frequency  0) 
and  the  equipment  (frequency  u) ) . K = & — co  ) /to  . 

B,  8 = damping  ratio  for  structure,  equipment,  respectively 
B = C/2ftM 

8 = c/2com 

j 

C,  c = viscous  damping  coefficients  for  structure  and  equipment, 

respectively . 

M,  m = masses  of  structure  and  equipment,  respectively 
y = m/M 
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In  order  to  apply  Equation  (2-1),  the  free-field  spectral  acceleration 

is  evaluated  at  the  shifted  frequency  (1  + and  at  the  average  damping 
ft  + R 

ratio  (-~2~ ) . Then  the  modification  factor  is  evaluated  which  depends  on 
the  detuning  £,  the  mass  ratio  y,  and  the  damping  ratios  3 and  B. 

In  the  event  that  the  mass  ratio  y « £2  + 43B  and  £ « 1,  corres- 
ponding to  light,  nearly  tuned  equipment.  Equation  (2-1)  becomes  equivalent 
to  the  conventional  floor  spectrum  approach  and  has  the  further  advantage  of 
being  applicable  when  the  only  available  information  on  ground  motion  is  a 
free-field  design  spectrum. 


SECTION  3 


MODAL  ANALYSIS  OF  STRUCTURE-EQUIPMENT  SYSTEM 


The  main  characteristic  of  the  interaction  between  the  structure  and  the 
tuned  equipment  is  in  the  equivalent  two-degree-of-f reedom  system  comprising 
the  equipment  and  the  particular  mode  of  the  structure  to  which  it  is  tuned. 
For  comparison  with  the  modal  analysis  of  the  structure  the  appropriate 
eauations  of  a two-degree-of-f reedom  system  are  given. 

The  equivalent  two  degree  of  freedom  system  is  shown  in  Fig.  2 in  which 

lower  case  letters  refer  to  the  equipment  and  upper  case  letters  to  the 

suitable  modal  properties  of  the  structure.  The  equivalent  or  effective 

ef  f 

ground  motion  is  denoted  by  u^  (t)  and  the  equations  of  motion  are: 


MU  + C(U  - ueff)  + K(U  - ueff)  = F(t) 
g g 

F(t)  = - m'u  = c(u  - U)  + k(u  - U) 


(3-1) 


where  F(t)  is  the  interaction  force  between  the  two  systems  and  where  m and 
M are  equipment  and  structure  modal  masses,  k and  K are  the  respective 
stiffnesses  and  c and  C are  the  damping  coefficients.  It  is  convenient 
to  introduce  the  variables 


oj  = / k/m 
(3  = c/2tom 


£2  = / K/M 
B = C/2f2M 


the  detuning  parameter  £ = (f2-w)/aj  and  the  mass  ratio  y = m/M  in  terms  of 

which  the  equations  become: 

" • 2 eff  2 eff 

U + 2BftU  + !1U  = 2Bfiu  + Q u + F/M 

g g 

F(t)/m  = -ii  = 20oj(u-U)  + u)^(u-U) 


(3-2) 


Applying  the  Laplace  transform  to  these  equations  and  eliminating  the 
structure  response  gives  a single  equation  for  the  transform  of  the  equipment 


1 


response  in  the  form 


C[(p2  + 2go)p  + oj2)  + lE^2&aL±J-2)  j . -eff  (3.3) 

p + 2B(lp  + fl  p + 2Bfip  + fl  ® 

The  equations  of  motion  of  the  N degree  of  freedom  structural 
system  take  the  form 

N ..  N , N N N 

l M..U.  + l C..U.  + T K .U.  = Y C . r .u  + 7 K..r.u  (3-4) 

j«l  13  j J-i  i3  3 j=l  13  3 j=l  13  3 8 j=l  13  3 § 

+ F(t)e. 

where  M_  is  the  mass  matrix,  C_  , the  damping  matrix  and  the  stiffness 

matrix.  The  vector  r is  a vector  of  influence  coefficents  introduced  to 

couple  the  actual  ground  motion,  u (t)  to  the  structure  and  e a vector 

g i 

which  is  zero  at  every  degree  of  freedom  except  that  one  denoted  by  the  index 
r,  to  which  the  equipment  is  attached  where  it  takes  unit  value.  F is  the 
interaction  force  between  the  equipment  and  the  structure. 


The  natural  frequencies  Q and  modes  $ are  given  by  the  equation 


l2 

n 


l I v” 


N N 

M $n  = 

j-i  13  J j-i 

Assuming  that  the  damping  is  sufficiently  small  that  it  does  not  introduce 
coupling  between  the  modes  the  equations  in  modal  coordinates  become 

\ + 2\flA  + - j1  <Fi/Mk 


where 


and 


N N 


"<. = ■ “A  -Ufi 

N N 

F,  = 7 C . r u + 7 K..r.u  + FeJ 

1 >1  ij  j 8 >1  lj  J g 1 


(3-5) 


(3-6) 
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r 


The  Laplace  transform  of  the  structure  response  U.(t)  is  given  by 

N J 

' W lf_ 

y <t>  $ f (p) 

N <Vi  j 

0 = 'j  — — 

1 k-1  ^(p2  + 2Bk^p  + (£) 


VP>  = (CuV  + + fei 


The  corresponding  equation  of  motion  for  the  equipment  displacement  u is 


-mu  = F(t)  = c(u  - U ) + K(u  - Ur) 


or  in  terms  of  Laplace  transforms 


2-  - o - 

-p  u = F(p)/m  = (28wp  + (i)  ) (u  - U^) 


From  Eq.  3.9  a relationship  between  u and  is  obtained  in  the  form 


(p2  + 28u)p  + cu2)u  - (2Bwp  + a)  )U 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


and  from  Eq.  3.7  this  can  be  written  as 

u (p2  + 20ajp  + a)2)  n N 

, » J.  *i  *r  [f‘i  + J.  <CitV  + 

= ( 2 3top  + a,2)  l ii 

k=i  ^(p  + 2Bknkp  + sy 

2- 

But  we  also  have  F = -mp  u . Using  this  to  eliminate  F in  the  above,  leads 
to  the  final  transformed  equation  for  the  equipment  response. 

2 ? N mp2  (2fki>p  + w2)$k<t>k 

U [( P + 280)p  + 0)  ) + l r ^ r?  ] 

k-1  \(p  + 2BknkP  + sy 

$k  y <j>k  y (C  r p + K r ) (3.11) 

, N r,f\  i iiV  i«,  f/ 

- (26w  + u>2)  l — 21 5 G 

k-1  \(P  + 2BkfikP  + (I2)  8 
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r k 2 r k 

We  note  that  I K.£<t>^  can  be  written  as  ftk  l ant^ 

1=1  1=1  N k 

the  assumption  of  small  damping  allows  the  representation  of  £ 0^$^  as 

r k — ^ 

2B  ft,  2,  MO..  Thus  the  solution  for  u for  the  multidegree  of  freedom 

tv  lv^_^  XX/  X 

system  takes  the  form  2 

2 2 2 N m<i>r  + **2) 

u [ ( p + 20u)p  + oj  ) + p l 2 2 

k=l  M^(p  + 2B^ftkP  + ftk> 

s i*\ ! + flk><2e"" + "Z) 

V i~  1 X/— 1 - 

= I 2 2 U 

k=l  Mk(pZ  + 2BkftkP  + ftk)  g 

In  performing  the  inversion  of  the  Laplace  transform  by  the  use  of 
residue  theory,  we  are  interested  in  the  zeros  of  the  term  in  brackets  on 
the  left  hand  side.  These  zeros  are  the  poles  of  the  transfer  function. 
Here  we  are  restricting  attention  to  the  case  where  the  equipment  frequency 


(3.12) 


is  close  to  a structural  frequency,  ft  , say.  This  is  indicated  in  Fig.  3. 

n 

In  this  figure  the  two  expressions  in  the  brackets  on  the  left  hand  side  of 
Eq.  3.12  have  been  plotted  separately.  These  plots  were  obtained  by  replacing 
p by  ift  and  then  drawing  the  graph  of  each  of  the  two  resulting  functions  in 
the  bracketed  expression.  For  simplicity  of  illustration  we  have  plotted  the 
figure  for  the  completely  undamped  case  (8  = = ...  = Bn  = 0).  The  negative 

of  the  first  function  plots  as  a simple  quadratic  in  ft,  becoming  zero  at  ft  = a), 
the  natural  frequency  of  the  equipment.  The  function  involving  the  summation 
plots  as  the  complicated  curve  which  goes  to  ± 00  at  ft  = ftk>  k = 1,  2....N, 
the  natural  frequencies  of  the  structure.  Two  such  curves  have  been  plotted, 
one  when  the  equipment  mass  is  small  and  one  when  it  is  not. 

The  values  of  ft  at  the  intersections  of  these  two  curves  give  the 
locations  of  the  zeros  of  the  bracketed  expression  on  the  left  hand  side  of 
Eq.  3.12,  which  are  the  poles  of  the  transfer  function  for  the  equipment 
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Location  of  Poles  of  Equipment  Response  Transfer  Function  for  Nearly  Tuned  NDOF  Structure 


response,  taking  into  account  equipment-structure  interaction.  It  is 

seen  that  in  the  case  of  small  equipment  mass  these  poles,  all  of  which  are 

simple,  occur  near  the  natural  frequencies  of  the  structure.  It  is  important 

to  note  that  two  closely  spaced  poles,  which  we  call  the  tuning  poles, 

are  located  near  the  equipment  frequency,  one  below  it  and  the  other  above 

it  as  shown  in  the  figure.  These  two  closely  spaced  poles  coalesce  into 

a double  pole  when  u)  = fi  and  m ->  0.  Thus  the  contribution  to  the  sum  of 

n 

the  residues  at  all  of  the  poles  is  dominated  by  those  which  are  associated 
with  the  two  tuning  poles.  In  computing  the  residues  at  the  tuning  poles,  it 
is  clear  that  the  contribution  from  the  summation  expression  is  dominated  by 
the  term  where  k = n because  then  the  denominator  of  that  term  is  nearly  zero. 
Hence  in  the  region  of  p = io),  Eq.  3.12  can  be  approximated  by 


- 2 2 2 m<I>n 
u[(p  + 2$wp  + a)  ) + p r 


M 


2go>p  + co . 

2 . _ „ . „2J 


n p + B p + ft 
r n nr  n 


(2B  ft  p + ft2)(28wp  + w2)  N N , n, n, . 

n n n'  r ' , r v $ $.Mjnrn/M  )u 
2 2 < Z Z r x iH  H n g 

p + 2B  S2  + Q i-1  £=1 

n n n 


(3.13) 


Comparing  this  expression  to  that  for  the  two-degree-of-f reedom  system, 
Eq.  3.3,  we  see  that  the  effective  mass  ratio  is 

Yeff  = m$n  /M 
r n 

and  the  equivalent  ground  motion  is 


(3.14) 


eff 


Cn  u 
r g 


where 


N N 

C"  = 0n  l l 0n  M.  r /M 

r r i-1  j-1  1 ^ J n 


(3.15) 
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The  solution  technique  used  here  will  be  to  obtain  the  contribution 
in  the  region  near  p = iw  by  considering  the  equivalent  two-degree-of-freedom 
system  defined  by  the  above  equations.  The  contribution  to  the  solution  from 


the  tuning  poles  requires  special  treatement,  and  this  will  be  done  in  the 
context  of  the  equivalent  two-degree-of-freedom  system.  The  contributions 
at  the  other  (N-l)  poles  is  straightforward  and  wil]  be  included  after  the 
two-degree-of-freedom  analysis  has  been  completed.  It  should  be  emphasized 
that  the  use  of  an  equivalent  two-degree-of-freedom  system  is  not  essential, 
but  is  only  conceptual  and  introduces  no  further  approximations  beyond  that 
made  in  passing  from  Eq.  3.12  to  Eq.  3.13. 


SECTION  4 

ANALYSIS  OF  TRANSFER  FUNCTION  FOR  TWO-DEGREE-OF-FREEDOM  NEARLY  TUNED  SYSTEM 


In  the  previous  section  the  parameters  which  appear  in  the  equivalent 

two-degree-of-freedom  system  have  been  derived.  Returning  to  the  Laplace 

transform  of  the  equation  of  motion  of  the  two-degree-of-freedom,  Eq.  3.3, 

to  simplify  the  subsequent  notation,  the  terms  B and  Q are  used  to  represent 

the  structural  parameters  B , of  the  tuned  mode  and  Y and  u should  be 

n n g 

interpreted  as  the  effective  mass  ratio  and  ground  motion  as  given  by 
Eqs.  3.14  and  3.15.  In  terms  of  these  parameters  the  transformed  equipment 


acceleration  u(p)  takes  the  form 


where 


u = [N(p)/D(p) ]u  (4. 1) 


N(p)  = (2(3u)p  + u>2)(2Bcd(1+£)p  + (l+£)2u)2)  (4.2) 


D(p)  = p4  + wp3(2B(l+0  + 2B(1+Q)+  o)2P2(2  + y (4>3) 

+ 2EJ  + C2  + 262B(l+fJ)  + oj3p(26(l+FJ2  + 2B(1H))  + W4(l+Q2 


where  E,  = 


(2  - u> 


is  the  detuning  parameter.  In  what  follows  attention 


will  be  concentrated  on  the  equipment  acceleration.  Completely  parallel 
results  can  be  easily  developed  for  the  equipment  displacement. 

The  nature  of  the  solution  depends  essentially  on  the  zeroes  of  the 
denominator  D(p) . Since  Y,  B,  B and  £ are  small  parameters,  the  roots 
of  D(p)  will  be  close  to  those  of  the  system  with  y»  3,  B and  5 taken  to 
be  zero;  namely 


p = + iw 
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To  determine  the  location  of  the  poles  of  D(p)  we  replace  p in  Eq.  4.3  by 

p = iw(l+«)  (4.4) 

where  6 is  a small  quantity.  Only  the  plus  sign  is  taken  since  the  roots 
will  appear  as  complex  conjugates.  In  terms  of  6 the  equation  D(p)  = 0 
takes  the  form 

<S4  + [4  - i[2B(l+y)  + 2B(l+£))]63 
+ [4  - y - 2£  - £2  - 4BB(l+£)  - 3i(2B(l+y)  + 2B(l+£))]<52 
+ [~2y  - 4£  - 2£2  - 8BB(l+£)  - i(28(2+3y  - 2£  - £2)  + 4B(l+£))]<$  (4.5) 

+ [-y  - 4$B(l+£)  - i2g(y  - 2£  - £2)]  = 0 


Solutions  for  this  equation  are  easy  to  obtain  when  B = 0,  B = 0 and y ^ 0, 
£ ^ 0 and  are 


2 2 2 4 ^ % 

6 = [l  + £ + -^  + -|  +(Y  + £2  + £3  + y£+  + + ] -1 

~ \ [£  +(£2  + Y)*0 


Also  if  y = 0,  £ = 0 and  8 ^ 0,  B ^ 0,  we  have 

6 = /l  - 8 2 - 1 + iB  , /l  - B2  - 1 + iB  £ iB,  iB  . 


Thus,  throughout  the  analysis  it  will  be  assumed  that  8,  B,  £ and  y5 
are  all  of  the  same  order,  say  e,  and  the  various  approximations  for  6 will 
be  based  on  assuming  6 of  order  £ <<  1 . When  the  parameters  are  not  of 

the  same  order  the  modifications  required  are  obvious. 


The  solution  of  Eq.  4.5  retaining  terms  of  order  c is 


6 = + — + i(fl  + ft  + >i) 

2 — 2 ' 2 — 2 


(4.6) 


\ 


i 
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where  here  and  throughout  the  remainder  of  the  analysis  the  upper  signs  are 
taken  together  to  give  one  root  and  the  lower  the  other.  The  quantities 
X and  Vi  are  given  by 

h ' * 

X = — [(tt2  + Y - (3-B)2)2  + 4£2(B-B)2|  + (£2  +y-(S-B)2)]  (4.7) 

fl  * > 

h *s 

U = — [|(C2  + Y - (3-B)2)2  + 4^2(6-B)2}2  ± (C2  +Y-(g-B)2)]  (4.8) 

It  is  easy  to  show  that  for  (3  ^ 0 and/or  B ^ 0 the  imaginary  part  is  always 
positive,  thus  leading  to  damped  oscillations.  In  view  of  the  large 
number  of  parameters  in  this  solution  there  are  many  special  cases  and  in 
the  following  sections  we  consider  some  of  these  of  particular  interest 
in  further  detail. 


4.1  UNDAMPED  TUNED  TWO-DEGREE-OF-FREEDOM  SYSTEM 

In  the  case  3=0,  B = 0 and  £ = 0 the  solution  of  Eq.  4-5  retain- 
3 

ing  terms  of  order  £ is 


(4.9) 


It  is  useful  in  this  case  to  retain  the  higher  order  terms  in  <5  since 
these  terms  are  necessary  in  deriving  certain  later  results  on  the 
floor  spectrum  method.  In  terms  of  the  transform  parameter  p the  roots 


p = iu)  + iw  — (1  ± •— — ) 

These  are  indicated  in  the  root  locus  diagram  Fig.  4,  with  the  corres- 
ponding complex  conjugate  roots.  It  is  clear  that  the  roots  remain  on 
the  imaginary  axis  with  the  small  spread  between  them  equal  to  Y"*  oj. 
These  will  lead  to  an  undamped  oscillating  solution. 


(4.10) 
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4.2  UNDAMPED  SLIGHTLY  DETUNED  TWO-DEGREE-OF-FREEDOM  SYSTEM 


In  this  case  the  roots  are 


H 


which  in  terms  of  p is 


<5  = f ± \ (Y  + K2) 


h 

p = i(D  (1  + -|)  + iw(y  + C2)/2 


(4.11) 


(4.12) 


(We  note  that  when  y -+  0 these  become  p = + ito,  + iS),  .)  These  roots 

and  the  corresponding  complex  conjugate  roots  are  indicated  in  Fig.  5. 

Again  an  undamped  oscillatory  solution  results  and  the  spread  between  the 

2 ^ 

closely  spaced  roots  is  now  given  by  to ( y + £ ) . 

4.3  DAMPED  TUNED  TWO-DEGREE-OF-FREEDOM  SYSTEM 

We  note  first  that  in  the  case  £,  = 0 there  is  the  possibility  of 
a double  root  of  the  equation  D(p)  = 0.  For  this  to  be  so,  certain  relation- 
ships must  exist  between  the  coefficients  of  the  various  powers  of  p in 
the  expression  for  D(p)  given  by  Eq.  18  when  £ = 0 . It  is  easy  to  show 
that  these  conditions  are: 

yS  = 0 (4.13) 

and 

2 7 9 

(4.14) 


y + 2BB  = B2  + (1  + y)2B2 


for  nonzero  y it  is  clear  from  Eq.  4.13  that  3 must  be  zero  and  from  Eq.  4.14 
2 

we  must  have  y = B 

*•  >. 

2 

The  solution  of  Eq.  4.5,  when  E,  = 0,  retaining  terms  of  order  e is: 

L 


6 = i ^T±  \ (Y  - (B-B)2)2 


(4.15) 
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Utilizing  this  as  a first  approximation  and  using  an  iterative  method 

3 

the  solution  of  Eq.  4.5  retaining  terms  of  order  £ is  obtained  in  the 


form: 


ft  n *1  2 ^ 

fi  = i - 2 — t - (3-B)  + (Y-(g-B)2)  (32+B2-y/2) 


+ i (26Y+(S+B)(Y-(e-B)  ))} 


(4.16) 


2 2 

When  the  term  Y - (3  _ B)  is  of  order  e it  dominates  in  the 
radical,  and  it  is  enough  to  use  Eq.  4.15  for  the  roots  6 . 


We  note  that  when  y -*•  0 this  gives  the  two  roots 


6 = i 3 and  6 = iB 


which  represents  the  floor  spectrum  solution  for  the  damped  system. 

2 

The  solution  (Eq.  4.15)  suggests  a double  root  when  Y = (3  - B)  , but  it 

has  already  been  shown  that  a double  root  will  only  occur  if  3 = 0;  thus 
2 3 

when  Y " (3  - B)  is  of  order  £ or  higher,  the  complete  expression  must 


be  used.  In  the  particular  case  y = (3  - B)  the  roots  are: 


= t 3 + b ± ± eV 


(4.17) 


When  3=0,  the  roots  become 


6 = i j±j  { (Y-B2)  + (Y-B2)Ss(B2-  + iB(Y-B2)} 


(4.18) 


For  fixed  B the  roots  go  from 


6 = i ^ j Y 2 f°r  Y » B2  to  6 = iB  - ~ , 0 for  Y "*■  0 


When  y = B , Eq.  4.18  indicates  a double  root  at 


<5  = iB/2 


However,  by  consideration  of  the  complete  Eq.  4.5  and  substitution  of  3=0, 
2 

Y = B it  is  easy  to  show  that  the  double  root  is  given  by: 
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(4.19) 


The  pattern  of  the  roots  in  the  p plane  for  3 = 0 is  shown  in  Fig.  6. 

2 

When  y ^ (6  - B)  the  nature  of  the  solution  depends  on  whether 
2 2 

Y > (3  - B)  or  y < (3  - B)  . For  the  first  the  roots  are  given  in  Eq.  4.15. 

For  the  latter  they  become,  to  lowest  order 

6 = i (3  + B + ((3-B)2  - Y)1S)/2  (4.20) 

The  roots  in  the  p plane  thus  have  the  same  imaginary  value  iw,  but  are 
equally  spaced  from  (3  + B)/2  on  a line  parallel  to  the  real  axis.  It  is 
important  to  note  that  both  roots  lie  in  the  left-hand  plane  for  all 
nonzero  values  of  3,  B and  Y so  that  the  transient  response  of  the  system 
always  decays. 

4.4  DAMPED,  SLIGHTLY  DETUNED  TWO-DEGREE-OF-FREEDOM  SYSTEM 

In  this  instance  the  large  number  of  parameters  makes  it  convenient 
to  illustrate  the  form  of  the  solution  by  considering  the  special  case 
when  6 = B.  Here 

2 h 

x = (y  + r) 

and  u = 0,  giving 

S - f + i(Y  + J2)*  ■*  <4-21> 

and 

£ 1 2 ^ 3+B 

p = iu>  (1  + | + j (y  + ? ) ) -w  — (4.22) 

These  roots  are  similar  to  those  shown  in  Fig.  5,  except  that  they  are 
shifted  to  the  left  by  <i>(-^^)  • 
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SECTION  5 


INVERSION  OF  TRANSFORM  SOLUTION  FOR  TWO-DEGREE-OF-FREEDOM  SYSTEM 


The  formal  inversion  of  the  transform  Eq . A. I is 


'd(t)  = 2¥T  J Vp)  ePt  dp 

r 

where  T is  a suitable  Bromwich  path.  If  u (p)  is  taken  to  be  1,  then  the 

8 

inversion  gives  directly  the  Green's  function  of  the  solution,  u„(t), 

G 

which  will  become  the  essential  ingredient  of  the  subsequent  analysis. 

The  complete  solution  for  the  acceleration  for  given  ground  motion  u (t) 
will  take  the  form 


(5.1) 


U ( t ) = U (t-T)  U (T)  dT 

G g 


(5.2) 


The  Green's  function  will  be  obtained  by  the  use  of  residue  theory,  there 
being  no  branch  cuts  in  the  p plane.  It  will  be  convenient  to  obtain  the 
inversion  of  the  transformed  Green's  function  for  the  general  case, 

Eq.  4.6  and  for  different  ranges  of  the  parameters  y,  6,  B and  £, 
corresponding  to  the  special  cases  discussed  in  detail  in  the  previous 


section. 


To  obtain  the  inversion  the  denominator  D(p)  is  written  in  the  form 


where 


D(p)  = (p  - PL)(P  - Pj)(p  - P2)(p  - P2) 


4 /,  , £ , A.  -6+B  . u, 

P1  = iu>  (1  + + j)  -0)  (— + j) 


and  p^  and  p^  are  the  complex  conjugates  of  p^  and  p^  • Evaluating  the 
residues  at  each  pole  and  collecting  complex  conjugate  terms  in  pairs 
lead  to  the  result,  correct  to  dominant  order. 
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••  / a.  \ u>  -oj(B+B)r/2  ( i , u A , . 

u (t)  = — — e l A sin  - ujt  cos  unn:  sin  w(i+  ^)t 

G A2  + p2  l 2 2 2 (5.3) 

- A cosh  w^t  sin  (jo~t  cos  w(l+f)t  - v sinh  -^uit  cos  cu|t  cos  w(l+^)t 

- p cosh  ^ ojt  sin  o)jt  sin  0)(l+|)t  | 

The  results  predicted  by  this  solution  are  explored  in  several  special  cases. 

5.1  Undamped,  tuned  system.  When  g = 0,  B = 0 and  K = 0,  the  solution 

3 

retaining  terms  to  the  order  e takes  the  form 

iir(t)  = ~ sin  u)t  cos  t|t  _ COs  out  sin  nt 

^ 2 y'* 

where  n = o>Y  /2  (5.4) 


5 . 2 Undamped,  slightly  detuned  two-degree-of-freedom  system . 

In  this  case  the  Green's  function  for  the  solution  takes  the  form, 

, 2 
retaining  here  terms  up  to  order  e , 

U(j(t)  — cos  w (l-t-|-)t  sin  r|t 

«T +y)"2 
2 1* 

where  now  H = (C  +Y)  2 0)/2. 


(5.5) 


5 . 3 Damped,  tuned  two-degree-of-freedom  system. 

For  the  situation  where  £ = 0 it  is  convenient  to  obtain  the  inversion 
of  the  transformed  Green's  function  for  three  different  ranges  of  the 
parameters  Y>  3,  B. 

2 

i)  For  y >(g  - B)  we  can  write  D(p)  in  the  form 

D(p)  = (p-iuH-Ejw)  (p+iuH-f yd)  (p-iu)  +e2u>)  (p+iort-t^oj) 

where 


£,  = 


g + B 


+ 2 (Y  - 


(3-B)2)’^ 


c 

1 


3 + B 
2 


(6-B)  )** 
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Evaluating  the  residues  at  each  pole  and  collecting  complex  conjugate 
terms  in  pairs  leads  to  the  result,  correct  to  uominant  order. 


uG(t)  = 


a)  e (B+B)u)t/2  cosu)t  sjn  ,/y.-(g-B) 2 cot/2 


(5.6) 


/RFb7 


This  function  represents  a damped  beat  type  solution,  the  beat  frequency 
being  w /y  - (3-B) r /2  which  is  much  smaller  than  the  tuning  frequency  u) 


ii)  For  y <(S_B)2  writing  D(p)  in  the  same  form  as  before,  evaluating 
the  residues  at  each  pole  and  collecting  terms  in  conjugate  pairs  leads 
to  a Green's  function  of  the  form: 


^ (i)  e (S+B)<i)t/2  cogwt  sinh  /(g-B)*  - y cot/ 2 

u„(t)  = - == 


/(B-B)  - 


(5.7) 


2 2 

Since  (B+B)  > (B-B)  - y for  non-zero  B,  B and  y>  the  term  exp[-(B+B)wt/2] 


dominates  the  term  sinh  / (B-B) * - y cut/ 2 . The  solution  can  be  interpreted  as 
overdamped  beats  by  analogy  with  the  concept  of  overdamped  vibrations.  For 
large  values  of  wt  the  solution  has  the  appearance  of  an  oscillation  of 
frequency  co  damped  by  an  exponential  with  factor 


- \ (B+B  - /(B-B)2  - Y))wt 


iii)  For  y = (B-B)  the  result  in  Eq.  4.15  predicts  a double 
pole.  However,  we  have  already  shown  that  if  B ^ 0,  a double  pole  will 
not  appear.  In  fact,  the  more  accurate  location  of  the  root  gives: 


p = io)  - — "-"B  u + i u)  S2Y^/2 


(5.8) 


Proceeding  in  the  same  manner  as  before  to  evaluate  the  residues  the 
following  result  is  obtained: 
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2 


uG(t)  = - 


C0  e 


-(B+B)u>t/2 


H h 

cosojt  sin  t3\  <i)t/ 2 


3V 


As  in  Eq.  5.6  this  represents  damped  beats.  The  beat  period  is  of  order 

_3  /2 

£ and  is  thus  very  long. 


iv)  When  3=0  and  y = B a genuine  double  root  appears.  The 
Green’s  function  solution  in  this  case  takes  the  form: 


uG(t)  = 


1 2 -Bait/ 2 

2 w t e cos  cot 


This  case  can  be  interpreted  as  critically  damped  beats. 


5.4  Damped,  slightly  detuned  two-degrees-of-freedom  system. 

In  the  special  case  of  a detuned  damped  system  with  3 = B the 
appropriate  Green's  functior  solution  is 

-oi(3+B)  t/2 


uG(t) 


(1) 


(Y+C2)4 


sin  r) t cos  ti)(l+£/2)t 


where  n = (y+£2)  * u)/2. 


(5.9) 


(5.10) 


(5.11) 
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SECTION  6 


APPLICATION  TO  EQUIPMENT  MOUNTING  DESIGN; 

UNDAMPED  TUNED  SYSTEMS 

The  results  obtained  in  the  previous  section,  for  the  response  of 
various  types  of  damped,  undamped,  tuned  or  untuned  systems  can  be  utilized 
in  the  design  of  equipment  or  equipment  mounting.  The  least  complicated 
forms  of  these  equations  are  those  for  the  undamped  tuned  system  and  it  is 
worthwhile  to  examine  these  in  detail  before  proceeding  with  more  general 
cases.  Many  of  the  basic  features  of  the  phenomena  are  more  readily 
apparent  for  this  situation.  The  methods  developed  for  the  damped  and 
tuned  or  detuned  cases  are  extensions  of  the  method  developed  for  this 
case . 

The  results  given  in  Eqs.  5-2  with  5-4  could  in  principle  be  used 
by  a designer  of  equipment  or  equipment  mounting,  if  a specified  ground 
acceleration  history  were  available  to  estimate  the  forces  which  would  be 
developed  in  the  equipment  or  its  mounting.  However,  such  information  is 
not  readily  available  to  a designer  and  the  computation  involved  in  these 
integrals  may  also  be  inconvenient  during  the  design  process.  It  would  be 
more  common  to  begin  with  a design  spectrum  which  may  be  specified  by  a 
code  or  determined  from  averaging  several  possible  inputs  as  for  example  in 
seismic  design.  Reference  14.  We  are  thus  interested  in  determining 
to  what  extent  the  results  in  Eqs.  5-2  with  5-4  can  be  used  to  provide 
estimates  of  maximum  acceleration  when  the  information  available  is  the 
response  spectrum  of  the  ground  motion  u^.  In  the  following  sections  a 
number  of  alternative  approaches  are  explored. 


H 
* | 

f 

* 

» 

I 

ft  j 
ft 


T 


6.1 


FLOOR  SPECTRUM  ANALYSIS 


When  the  Green's  function  obtained  for  the  undamped  tuned  system 
Eq . 5-4  is  substituted  in  Eq . 5-2  the  response  is  given  by 

U(t)  = | iig  (t ) | | Y2  cos  n(t-T)  sin  a)(t-T) 

Y 0 

- sin  n(t-T)  cos  w(t-T)|  dt 


with  n = ojY  /2. 

For  small  values  of  nt  this  reduces  to 

’u(t)  = io  | u^(t)  | 2 sin  ~ w(t-T)  cos  u>(t-T)l  dt 

0 

Note  that  this  result  is  independent  of  y and  could  be  obtained  directly 
by  means  of  the  floor  spectrum  analysis  method  whereby  the  input  to  the 
structure  is  used  to  compute  the  base  motion  at  the  equipment  assuming  the 
equipment  to  be  absent  and  the  equipment  motion  is  calculated  with  this 
motion  as  input.  It  is  clear  that  this  approximation  is  valid  for  small 
values  of  nt,  i.e.  ojY  t/2  l.  The  essential  characteristic  of  this 
result  is  that  it  neglects  interaction  between  the  structure  and  the 
equipment . 

One  other  feature  of  the  floor  spectrum  analysis  is  that  it  can 
be  obtained  from  the  basic  equations  by  setting  y = 0.  A double  pole  will 
then  appear  in  the  tuned  case  (for  an  untuned  system  only  simple  poles 
occur).  This  double  pole  leads  to  terms  in  t cos  0Jt  and  t sin  U)t  on 
inversion.  Thus  the  floor  spectrum  analysis  cannot  be  used  to  determine 
maximum  displacement  or  acceleration  for  undamped  tuned  systems  since  it 
yields  responses  which  grow  without  limit.  It  follows  that  although  the 
floor  spectrum  analysis  is  a valuable  method  for  untuned  systems  it  has  no 
meaning  for  undamped  tuned  systems. 


(6-1) 


(6-2) 
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MODIFIED  GROUND  MOTION  SPECTRA 


When  nt  is  not  much  smaller  than  unity  the  term  in  Eq . 6-1 

j, 

multiplied  by  y2  is  negligible  in  comparison  to  the  other  and  it  becomes 


u(t)  = - a)  J 'u^(x)  sin  n(t-0  cos  w(t-r)dx 
0 

Expanding  the  term  cos  ca(t-x)  in  Eq . 6-3  allows  it  to  be  written  in  the 


(6-3) 


ii(t)  = - to  cos  (tot  — 4>) 


) | (f  iaf 

* n 


(t)  cos  cot!  sin  q(t-x)dx) 


+ (j  [up(T)  sin  <oT  ] sin  n(t-l)dT) 


(6-4) 


where 


f u t ( x ) sin  ojx]  sin  n(t-x)dx 


[u  (x)  cos  cox]  sin  r)(t-x)dx 


(6-5) 


The  terms 


— [u  (x)  cos  to"!)  sin  H(t-x)dx 
* ) 8 
0 


[u  (x)  sin  ojt]  sin  q(t-x)dx 


can  be  interpreted  as  the  response  of  a singl e-degree-of-f reedom  system 

with  frequency  q to  the  modified  ground  input  accelerations  u (t)cos  tot 

8 

and  u (t)  sin  tit.  We  aiso  note  that  n <<  to.  The  term  cos  (ait-if)  is  a 
g 

rapidly  oscillating  function  and  achieves  its  maximum  many  times.  The 
integrals  are  slowly  oscillating  functions  and  represent  a slowly  varying 
envelope  of  the  more  rapidly  oscillating  term.  The  maximum  of  the  product 
is  accordingly  very  nearly  the  maximum  of  the  envelope. 


Thus,  one  way  to  estimate  the  motion  of  the  equipment  is  to 


construct  spectra  for  the  modified  ground  accelerations.  From  this  an 
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estimate  of  the  maximum  acceleration  is 


u = |s^(n)2  + s^(n)2  l 

max  | A A f 


1/2 


(6-8) 


A similar  expression  can  be  obtained  for  the  displacement.  In  the  above, 
the  terms  S^(.'l)  and  S^(n)  are  undamped  acceleration  response  spectra  for 

A A 

the  modified  ground  motions  u (t)  cos  U)t  and  ii  (t)  sin  cot,  respectively. 

8 8 

In  principle  then  it  is  possible  to  develop  a design  technique  if  the  time 
history  ii^(t)  *s  available  by  constructing  the  low-frequency  response 
spectra  for  the  modified  ground  motions.  In  certain  cases  the  given 
information  may  be  only  the  response  spectrum  of  u^  and  not  the  ground 
motion  itself.  As  far  as  can  be  seen  at  the  moment  there  is  no  way  to 
compute  the  spectra  for  the  modified  ground  motion  if  the  only  information 
available  is  the  spectrum  of  the  actual  ground  motion. 


6.3 


AMPLIFIED  GROUND  MOTION  SPECTRUM 


In  view  of  the  remarks  in  the  previous  section,  we  now  develop  an 
alternative  approach  in  which  the  term  sin  D(t-x)  is  expanded,  leading  to 

2 


u(t)  r sin  (nt  - 0) 

y"5 


{<f\ 

' 0 


(x)  cos  r)T  cos  co(t-T)dT) 


+ ( 


2 1 ^ 


ii  (x)  sin  ryu  cos  w(t-x)dx) 

g 


(6-9) 


where 


0 
r t 


0 = tan 


-1  0 


ii  (x)  sin  r|T  cos  u)(t-x)dT 

g 


(6-10) 


ii  (t)  cos  HT  cos  oj(t-x)dx 
g 


We  are  interested  in  situations  where  the  ground  motion  has  a prescribed 
finite  duration  and  fpj  those  frequencies  ui  where  the  maximum  response  of  a 
single-degree-of-f reedom  oscillator,  i.e . the  response  spectra,  is  achieved 
late  in  or  after  the  termination  of  the  ground  motion.  These  frequencies 
correspond  to  peaks  in  the  response  spectrum  of  ground  motion  of  the 
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earthquake  type.  Design  spectra,  reflecting  the  probabilistic  motion  of 
the  input,  correspond  closely  to  the  peaks  of  actual  spectra  and  thus 
presuppose  late  occuring  maxima.  In  blast-type  ground  motion  which  is  of 
short  duration  it  is  likely  that  the  maximum  values  of  equipment  response 
will  occur  at  times  larger  than  the  duration  of  the  ground  motion. 


Thus,  for  values  of  f|t^  <<  1,  where  t^  is  the  duration  of  the 
ground  motion,  the  first  integral  in  Eq.  6-9  can  be  approximated  by 


ft 

’u  (x)  cos  io(t-x)dx 

g 

0 

and  the  second  neglected  since  sin  nt  will  be  bounded  by  r)t^  « 1.  For 
Ht1  <<  1,  then,  we  have 


**  Cl) 

u(t)  = - sin  nt  I ii  (t)  cos  (a(t-x)dx 

v'5  J g 

y 0 

The  term  in  the  integral  is  a function  oscillating  with  frequency  co 
which  is  high  compared  to  n and  a maximum  of  that  term  will  nearly 
coincide  with  the  maximum  of  sin  X)t . An  estimate  of  the  maximum  value 
of  ii(t)  is 


max 


u (x)  cos  co(t-x)dx 

o 


If  the  displacement,  velocity  and  acceleration  response  spectra  as 
functions  of  frequency  w and  damping  parameter  6 are  denoted  by  S^(u),  B)  , 
Sv(ia,  B)  and  SA(w,  B)  respectively,  we  then  recognize  that 


ft 

max 

11  (x)  cos  co(t-x)dx 

1 g 

0 


is  the  undamped  velocity  response  spectrum  Sv(u),  0)  for  a single-degree- 
of-freedom  system  with  frequency  u.  Thus,  we  have  as  an  estimate  of  u 

max 


to 


S„(u,  0) 


the  expression 


Since  = Sy/cu  and  = u>Sy,  then 

S_(a>,  a) 

|u|  = -±~r  - 

max  ‘2 

Y 

Sa(oj,  0) 

max  h 

Y 

It  follows  that  if  an  engineering  designer  is  given  only  the  response 
spectrum  of  the  ground  motion,  the  maximum  displacement  and  force  in  the 
equipment  can  be  estimated  by  using  these  spectra  amplified  by  the  factor 
Y 2-  These  remarks  refer  of  course  to  the  equivalent  two-degree-of-f reedom 


system.  The  results  for  the  general  system  are  obtained  by  utilizing  the 
factors  in  Kqs.  3-14  and  3-15. 


SECTION  7 


APPLICATION  TO  EQUIPMENT  MOUNTING  DESIGN; 

DAMPED,  TUNED  AND  SLIGHTLY  UNTUNED  SYSTEMS 

In  the  previous  section,  a number  of  procedures  were  developed  for 
undamped  tuned  systems  which  could  provide  information  to  a designer  of 
equipment  or  equipment  mounting.  These  procedures  were  motivated  by  the 
fact  that  in  many  cases  a specified  time  history  of  the  ground  motion 
applied  to  the  structure  would  not  be  available  except  in  the  restricted 
form  of  a design  spectrum.  These  results  enable  the  designer  to  utilize  the 
design  spectrum  for  the  structure  to  estimate  directly  maximum  values  of  the 
equipment  acceleration  and  displacement. 

Three  different  approaches  to  the  design  problem  were  developed. 
These  were  the  floor  spectrum  method,  the  modified  ground  motion  spectrum 
method  and  the  amplified  ground  motion  spectrum  method.  It  was  shown  that 
for  tuned,  undamped  systems  the  floor  spectrum  method  was  not  a valid 
technique.  The  modified  ground  motion  spectrum  method  was  valid,  but 
inconvenient  and  it  is  clear  that  the  amplified  ground  motion  spectrum  was 
the  most  convenient  to  use  in  estimating  the  response  of  light  undamped 
equipment. 


For  damped  tuned  systems  it  is  seen  from  Eq.  4-3  that  if  7 <<  1, 
then  y is  negligible  if  4$B  >>  y.  This  means  that  for  such  cases  the 
possibility  of  significant  interaction  between  equipment  and  structure  can 
be  ignored  and  the  floor  spectrum  method  used  to  determine  the  response  of 
the  equipment.  However,  the  floor  spectrum  method  requires  the 
computation  of  time  histories,  and  thus,  if  a design  spectrum  is  the  only 
given  information,  its  use  may  not  be  the  most  convenient. 

Since  the  results  of  the  previous  section  showed  that  the 
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amplified  response  spectrum  method  was  the  most  convenient  to  use,  similar 
approximations  will  be  developed  for  damped  and  slightly  detuned  systems. 


UNDAMPED,  SLIGHTLY  DETUNED  SYSTEMS 


The  acceleration  response  u(t)  to  imposed  ground  acceleration 

u (t)  is  obtained  from  Eqs.  5-2  with  5-5  in  the  form 
g 


ii(t)  = - 


— r ii  (t)  sin  n(t-T)  cos  u)(l  + J)  (t-t)dT 

(£  +Y)"5  i 8 2 


(7.1) 


where 


n = (S2+y)'V2 


Expanding  the  term  sin  n(t-T)  as  was  done  in  Eq.  6-9  of  the 
previous  section  and  neglecting  the  analogous  terms  an  estimate  of  the  maximum 


acceleration  is  obtained  as 


SA  (u)(l+£/2),0) 

(yH  ) 2 


(7.2) 


This  result  is  still  valid  if  y <<  £ providing  £,  <<  1.  For  such  cases  the 
floor  spectrum  method  is  applicable,  but  of  course  could  not  be  used  if  the  only 
information  on  the  ground  motion  is  a design  spectrum.  The  above  result  is 
clearly  more  convenient  and  equally  valid.  The  beating  phenomenon  which  is 
the  physical  basis  of  the  result  will  appear  in  the  floor  spectrum  solution 
in  the  slightly  detuned  case,  the  beating  being  between  the  two  closely 
spaced  frequencies  U)  and  a)(l+£). 


DAMPED,  TUNED  SYSTEMS 


The  results  obtained  in  Section  3 for  the  transfer  functions  for 
damped  tuned  systems,  indicate  that  four  sets  of  the  parameters  y,  3 and  B 
have  to  be  identified. 
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Case  1:  y > (3-B)‘ 


The  acceleration  response  ii(t)  to  imposed  ground  acceleration 


ii  (t)  is  given  by 
8 


r t 


u(t) 


where 


(Y-(R-B)2)?5 


n = | (y-(6-b)2)^ 


ii  (t)  e t)/2  cos  U)(t-'t)  sin  n(t-x)dx 

g 


(7.3) 


Expanding  the  term  sin  n(t-x)  leads  to 


ii(t)  = — cos(nt-O) 

(Y-(f3-B)V2 


r t 


ii  (t)  e 
g 


-(3+B)w(t-x)/2 


0 


cos  rix  dT ) + ( 


ft  ••  , , - ($+B)  o)(t-T)  /2  . . . . . 

u (l)  e cos  0)(t-T)  sin  n't  di} 

8 


cos  0)(t-T) 
2/  *5 


where 


0 = tan 


-1  J 0 


ii  (T)  cos  nx  e-(8+B)u(t-T)/2cos  a»(t-T)d/ 
g 


ii  (t)  sin  px  e 

O 


-(g+B)ai(t-T)/2 


cos  co(t-x)dX 


0 


7.4) 


(7.3) 


We  are  interested  in  situations  where  the  ground  motion  has  a prescribed 
finite  duration  and  frequencies  for  which  the  maximum  response  of  a single- 
degree-of-f reedom  oscillator  is  achieved  late  in  or  after  the  termination  of 
the  ground  motion.  Thus,  for  values  of  r|t^  <<  1 where  t^  is  the  duration  of 
the  ground  motion  the  first  integral  in  Eq.  7-4  can  be  approximated  by 

(t  ..  , , - (B+B)co(t-T)/2  . , 

u (x)  e cos u)( t-  t) d t 

g 

0 

and  the  second  neglected  since  sin  nx  will  be  bounded  by  nt^  <<  1 and  ii^  = 0 
for  t > t^.  Thus  we  take 


, w sin  nt 

u(t)  '2'Y. 

(y-(3-b)  ) 2 


ft 


••  ,rs  -(6+B)w(t-x)/2 

ug(r)  e cos  to(t-T)dt 


(7.6) 
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When  the  parameters  y*  8 and  B are  small  we  may  interpret  this  result  in 
the  following  way: 

For  t > tL  the  above  expression  can  be  written  in  the  form; 


^(t)  = _ -(e+B)wt/2  _ 

2n  R 


cos(tot-'l') 


where 


and 


R = 

(A^ 

2,JS 
+ A ) 

l 

ft  - 

2 

A = 

1 

u ( t ) e 

1 

J 

0 

ft. 

g 

A = 

1 

u (t)  e 

2 

J 

0 

g 

V = 

tan 

< 

< 

with 


sin  tot  dt 


In  the  above,  R and  'V  are  Constants  independent  of  t for  t > t , and 
R cos (tot-f)  is  a rapidly  varying  function  of  time.  The  term 

^2si»  nt  e-<B«»»t/2 

is  a slowly  varying  envelope  curve  which  attains  its  maximum  value  at  a time 
t*  given  by 


tan  nt*  = 2n/(8+B)to 

The  value  of  sin  nt  at  which  the  envelope  achieves  its  maximum  is 
sin  nt*  D 

It  follows  that 


(7.7) 


(n'+(6+B)  Vm)"5 


(7.8) 


|lJ,max  = = 


to 


(y-(8-B)?) 


L I sin  nt* 


ft* 

'1  “ 


,T)  -(B+B)io(t*-x ) /2  , x 

,(T7  e cos  0)(t*-T)dT 


(7.9) 
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Note  that  the  term 


C ..  . . -(&4-B)u>(  t-T)  /2  , . 

u (t)  e cos  u>(t-T)dT 

g max 


8+B. 


is,  to  the  order  of  8 and  B,  the  velocity  response  spectrum  S y(w,  — — ) for 
a damped  single-degree-of-freedom  oscillator  with  damping  factor  (8+B)/2  and  is  a 
bound  for  the  integral  in  Eq.  7-9.  Thus  an  estimate  for  the  maximum  equip- 
ment acceleration  is 

|u|  . "I-8.1?  s g|S) 

Utilizing  the  value  of  sin  r)t*  from  Eq.  7-8,  the  final  estimate  is 


„ / B+B. 
w Sv^a)’  ~~2~ * 


' max 


(Y+46B) 

For  the  lightly  damped  systems  considered  here 


(7. 10) 


a)  sv  = SA  = 0)  SD 


so  that  the  result  can  be  written  in  the  alternative  forms, 


SA(“-  ifW48B>1/2 


(7.11) 


or 

V“-¥>/(2+''bb|1/2  <? 

It  follows  that  if  a designer  is  given  only  response  spectra,  at  various 
damping  values,  of  the  ground  motion  applied  to  the  structure,  the  maximum 
displacement  and  force  in  the  equipment  can  be  estimated  by  using  the 
appropriate  damped  spectra  for  a damping  factor  equal  to  the  average  of 
those  in  the  structure  and  equipment,  amplified  by  the  factor  (y+48B)  2. 

It  is  to  be  noted  that  if  8+B  is  fixed,  the  maximum  value  of  43B  is  achieved 
when  8 = B,  yielding  the  smallest  value  of  the  amplification  factor.  Thus 
if  the  total  damping  is  fixed  the  optimal  choice  is  to  have  it  shared 
equally  by  equipment  and  structure. 
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f 


Case  2:  y < (3-B)‘ 


The  solution  in  this  case  for  u(t)  in  terms  of  specified  ground 


motion  ii  (t)  takes  the  form 
8 


ii(  t)  = - 


((6-B)  -yr  l 8 


(T)  e“(6+B)  <i)(t'T)/2cos  w(t-r)  sinh  n(t-x)dx 


where 


(7.13) 


n = «B-B)2  - y )*  w/2 


This  can  be  written  in  a form  analogous  to  that  in  Eq.  7-4  with  the  envelope 

c -(6+B)  a)t/2  . , , , , 

now  in  the  form  e sinh  nt.  When  the  envelope  is  analyzed  as 

before  for  its  maximum  value  it  is  found  that  the  time  t = t*  is  such  that 


sinh  nt*  = n /(n2  + (&^b> 2 w2/4)is 


Following  the  arguments  used  to  obtain  the  previous  results  it  is 

k 

found  that  as  before  the  amplication  factor  is  (y+4BB)  . 


Case  3: y = (B-B)' 


The  solution  in  this  case  takes  the  form 

ft 


to 


u(t)  = “ 1 ft 


Ug (t)  e T)/2  cos  sin  ^(t-Tjdx 


(7.14) 


0 


where  now 


bV5/; 


n = (i)B"Tz/2 


In  this  case  the  envelope  takes  the  form 


sin  nte 


■ (B+B)u)t/2 


and  the  time  t*  is  given  by 


tan  nt*  = 2n/(3+B)w 
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Using  this  and  the  appropriate  value  for  n gives  as  the  amplification  factor 


the  term 


((0+B)2+By>'5 

3 

The  term  By  is  of  order  £ and  may  be  neglected  in  comparison  to 

2 -1 
(B+B)  . Thus  the  amplification  factor  is  simply  (&fB)  , but  note  that 

2 —1 

since  y = (8-B)  that  (6  +B)  = (y*-48B)  2.  Thus,  the  same  result  as  in  the  two 

previous  cases  applies  here. 


Case  4:  y = B , 8 = 0 


This  is  the  situation  in  which  a true  double  root  appears  and  the 
solution  has  the  form  of  a damped  floor  spectrum  result. 


u(t)  = - 


2 J g 
0 


(t)  e'Ba,^'T^2cos  w(t-x)  • (t-T)dx 


(7.15) 


The  envelope  is  wt  e ®wt^2  tjle  maxfmum  Qf  which  is  reached  at  t*  given  by 

tot*  = 2/B.  This  leads  to  an  amplification  factor  1/B.  We  note,  however, 
2 

that  since  y = B and  8=0,  the  amplification  factor  is  again 


(Y+48B) 


Thus  the  Eqs.  7-11  and  7-12  obtained  for  Case  1 are  in  fact  correct  for  all 
combinations  of  y,  8 and  B,  which  is  a surprising  result  when  one  considers 
the  differences  in  form  of  the  Green’s  functions  for  each  case. 


DAMPED,  SLIGHTLY  DETUNED  CASE 


The  addition  of  slight  detuning  considerably  modified  the  form  of 
the  response  which  now  becomes 
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U(t) 


to  (l  -<o(frhB)(t-T)/2 


2 7~ 

(A  +p  ) 


u^(t)  I - A sinh  ^ <o(t-x)  cos  w-(t-T) 
sin  io(l+|)(t-X) 


u A r 

+ A cosh  y w(t-x)  sin  uo^  (t-x)  cos  W(l+J)(t-T) 


+ p sinh  ^ co(t-x)  cos  10^  (t-x)  cos  to(l-t-^)  (t-x) 


(7.16) 


+ p cosh  ~ co(t-x)  sin  to--  (t-x)  sin  to(l+|)  (t-T)  } dx , 


where  A and  p are  defined  in  terms  of  £,  y,  3,  B in  Eqs.  3-7  and  3-8.  To 
simplify  the  algebraic  manipulations  needed  for  this  considerably  more 
complicated  expression  attention  is  focused  on  a single  case  which  will  be 
illustrative  of  the  result.  The  case  selected  is  that  of  optimal  use  of 
damping.  That  is,  if  a fixed  total  amount  of  damping  3 + B is  specified  then 
the  best  selection  of  this  damping  is  the  case  3 = B.  Then  Eq.  7-16  takes 
the  form 


u(t)  = - 


-co(3+B)(t-x)/2  ..  , . A . . £ ...  . . 

u (x)  sin  air  (t-x)  cos  (p(H-^)(t-x)dx 

g 2 2 


(7.17) 


where 


* ■ «2+y)1/2 


Expanding  the  sin  term  and  recalling  that  for  t >>  t^,  the  duration 

of  the  ground  motion,  the  term  sin  ~ cot,  can  be  neglected,  we  have 

u(t)  = - j e (3+B)u  I Acot/2  R cos(tot-t|c) 


where  here 


R = 


2 2 1/2 
(Aj+Aj)  ' 


and  ’I)  = tan  k^/k^ 


with 
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SECTION  8 


COMPLETE  SOLUTION  INCLUDING  OTHER  POLES 


In  the  previous  sections  we  have  determined  the  contribution  to 
the  response  from  the  tuning  poles  of  the  equipment-structure  system.  This 
is  of  course  the  dominant  part  of  the  response  in  the  case  of  light 
equipment  mass,  but  it  is  easy  to  include  the  contributions  from  the  other 
poles.  To  do  this  recall  that  the  non-tuning  poles  of  Eq.  3-12  are  close  to 
their  location  for  the  structure  alone,  as  indicated  in  Fig.  3.  For  the 
m*”*1  non-tuned  mode  the  poles  are 


p = - B fi  ± ifi 
mm  m 


Evaluating  the  residues  and  dropping  the  negligible  terms,  which  are  those 

k2 

multiplied  by  the  small  modal  mass  ratio,  m<J>r  /M^  <<  1,  k = 1 to  N,  we 
obtain,  to  dominant  order,  the  contribution  from  the  m*1*1  poles  as 


. - (fi  / to) 
m / 


-B  fi  t 

- fi  e m m sin  fit  ,m  4 n 
2 m m 


and  contributions  of  the  same  order  from  the  tuning  poles  as 


N C_ 

l 

m=l  l-(to/fi  )2 
m#n  m 


-Btot 

we  sm  tot 


where  C^  is  defined  in  Eq . 3-15.  Thus,  the  complete  solution  for  the 
response  of  the  equipment  takes  the  form 


ft  INC  , 

u(t)  - u (T)/  l -r 7l 

L 8 Jm-l  l-(fi  /to)  ' 

° (m*n 


/ "Bmnm  (t_T)  \ 

~z  ( fi  e sin  fi  (t— t ) ) 

, \ m m / 


+ I ~ ( toe  T^sin  to(t-x)  | + ii_(t-T)  1 dT 

m*l  l-(to/fi  )2  l / G ) 

m^n  m ' ' 


(8-2a) 


(8-2b) 


where  Ug(t)  is  the  appropriate  form  of  contribution  from  the  tuning  poles 
given  in  its  various  forms  in  Section  5. 

In  utilizing  this  complete  solution  to  develop  the  spectral 
response  in  each  mode  for  design  purposes,  it  is  important  to  note  that  the 
two  parts  of  the  solution  in  Eq.  8-3  have  an  entirely  different  character.  The 
contributions  from  the  non-tuning  poles  and  the  nondominant  contributions  from 
the  tuning  poles  are  conventional  and  would  attain  their  peaks  during  the  ground 
excitation  or  shortly  thereafter,  while  the  dominant  response  from  the  tuning 
poles  as  indicated  in  Sections  6 and  7 is  controlled  by  the  energy  transfer  from 
the  structure  to  the  equipment  through  beating,  which  takes  a relatively  long 
time.  The  maximum  in  the  latter  case  will  be  achieved  considerably  later  than 
the  former.  It  does  not  make  sense  to  add  these  in  the  conventional  way  such 
as  square  root  of  the  sum  of  squares,  or  by  a similar  rule.  In  fact,  they 
should  not  be  added  at  all,  but  treated  as  separate  maxima.  The  maximum  res- 
ponse from  the  nondominant  contributions  can  be  estimated  by  the  conventional 
method  of  square  root  of  the  sum  of  squares. 

Accordingly,  the  estimate  of  the  maximum  acceleration  has  two  parts, 
an  early  peak  given  by 


1 max 


« r 

1 h 

pi  Li- 


I m; 

. m^n 


N 

l 


(ft  /w) 

m 


O SA^  >B  ) 

2 A m*  m 


]■ 


(8-4) 


2)h 


m=l  1- (w/ftj 
m^n 


Sa(w,B) 


2 A 


m 


and  the  other,  a later  peak,  from  the  dominant  contribtions  of  the  tuning  poles 
given  by 


„n 


max 


/f2+Yeff+4BB^yi 


/ F \ 

^u^),  -2-°) 


(8-5) 


50 


where  £ = (12^-oj) /(a,  and  y is  given  by  Eq.  3-14.  For  light  equipment  mass 
and  lightly  damped  closely  tuned  systems  the  second  peak  will  be  the  most 
important . 


Although  not  of  immediate  interest  in  this  report,  it  is  also 
possible  to  utilize  the  methods  developed  to  obtain  estimates  of  the  peak 
response  for  systems  which  are  grossly  detuned,  i.e.,  where  the  equipment 
frequency  is  spaced  between  and  well  away  from  all  of  the  structural  frequencies. 


To  do  this  note  that  for  light  mass  the  structure  poles  are  only 
slightly  shifted  from  their  location  for  the  structure  alone,  namely 

p = - B 12  ± i!2m 

m m 


and  additional  poles  at 


p = - goo  + iw 


due  to  the  equipment  are  included  as  shown  in  Fig.  7 (illustrating  the 
undamped  case).  The  residues  at  the  structure  poles  are  as 
before  with  the  contributions  from  each  m = 1 to  N poles  given  by  Eq.  8-2a. 
The  residues  at  the  equipment  poles,  Eq.  8-7,  provide  a contribution  to  the 
Green's  function  in  the  form,  similar  to  Eq.  8-2b 


N Cm 

l ~ ~ 

m=l  1- (eo/Q  )2 

_ m _ 


-Scot  . 

to  e sin  cat 


The  derivation  is  completely  standard  and  is  similar  to  the  terms  from  the 
structure  poles.  The  complete  response  for  the  equipment  in  the  grossly 
untuned  case  is  thus  given  by 


(8-6) 


(8-8) 


fi(t)  = u (T)J  l -r 12  e m m 

I 8 |m=l  1—  (12  /u>)2  m 


-B  12  (t-x) 


Jin  12  (t-T) 
m 


/ N Cm  \ 

l 1 ) 

\m=  1 1-  (cu/12  )7 


-3td  (t-x)  . . J , 

0)  e sin  (J'(t-T)}  dT 


(8-9) 
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and  the  appropriate  estimate  of  the  peak  response  is  the  conventional  one 
using  square  root  of  the  sum  of  squares 


max 


S.(ft  , B ) 
A m m 


■] 


2 


C r Cr 

- 

2 

l r 

S (w,B) 

1 m=  1 l-;<i)/fi  )? 

A 

L m 

* 

, 

This  result  can  be  used  as  an  alternative  to  time  history  or  modal  analysis 
of  the  composite  N + 1 system.  It  is  also  an  alternative  to  the  standard  floor 
spectrum  analysis  which  requires  computation  of  the  time  history  of  the 
structure  alone.  This  is  next  applied  as  input  to  the  equipment  and  then  the 
time  history  of  the  equipment  is  computed.  Note  that  Eq . 8-9  is  completely 
independent  of  the  modal  mass  ratios.  In  fact,  it  represents  the  general 
closed  form  solution  of  the  floor-spectrum  method,  and  its  interpretation 
directly  provides  the  simple  estimate,  Eq.  8-10.  Indeed,  the  preceding 
analysis  which  led  to  Eq . 8-9  is  the  mathematical  justification  for  the  use 
of  the  floor  spectrum  method  for  the  grossly  detuned  system.  All  of  the 
information  needed  for  Eq . 8-10  is  available  from  the  building  design,  the 
equipment  frequency,  and  the  design  spectrum;  it  should  thus  be  very 
convenient  for  the  designer  to  use  in  practical  design  applications. 


(8-10) 


It  is  worthwhile  noting  that  the  methods  developed  for  dealing  with 
the  tuned  poles  can  be  used  to  determine  the  response  of  systems  with  closely 
spaced  modes  even  if  no  equipment  is  included.  Of  course,  in  this  case  the 
approximations  used  here  based  on  the  small  mass  ratio  could  not  be  used,  but 
the  treatment  of  the  envelope  of  the  beating  response  would  be  entirely 
similar.  It  is  not  to  be  expected,  however,  that,  in  the  case  of  closely 
spaced  modes  the  maximum  response  would  be  much  different  from  that  of  the 
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other  modes,  except  that  it  will  occur  at  a much  later  time.  It  is  the 


< 


light  equipment  mass  that  produces  the  large  amplification,  and  the  dominance 
of  the  result  in  Eq.  8-5  for  the  late  peak  over  that  in  Eq.  8-4  for  the  early 
peak. 

This  is  the  reason  that  the  peculiar  ad  hoc  approach  used  by 
Penzien  and  Chopra  in  Reference  11  led  to  a good  result.  In  their  approach, 
they  model  an  N-degree-of-f reedom  system  with  a light  appendage  by  consider- 
ing a set  of  N-two-degree-of-freedom  systems  in  which  one  component  is  one  of 
the  N modes  of  the  structure  without  an  appendage  and  the  other  is  always  the 
appendage  and  numerically  solve  the  N sets  of  coupled  differential  equations. 
They  then  add  together  the  peaks  from  each  of  the  two-degree-of-f reedom  systems 
by  the  square  root  of  the  sum  of  the  squares  to  obtain  the  maximum  response  of 
the  appendage.  That  the  result  is  reasonable  arises  from  the  fact  that  in  a 
tuned  system  the  contribution  from  the  tuning  poles,  and  thus  in  their  case 
from  the  particular  two-degree-of-f reedom  system  which  is  tuned,  dominates  the 
rest.  For  this  reason,  their  result  is  fortuitous.  Were  this  approach  used 
for  an  untuned  system,  it  could  produce  erroneous  results. 
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NOTATION 


V A2 


C 

Cij 

cm 

J 

D(  ) 

6i 
F(  ) 

j > k * ^ > 

etc . 
k 
K 


Kii 


M 

M.  . 
ij 

M 

n 

N 

N(  ) 
P 

q 


constants 

equipment  damping  coefficient 

model  damping  coefficient 

structure  damping  matrix 

modal  participation  factor 

transform  transfer  function  denominator 

equipment  attachment  vector 

interaction  force  between  equipment  and  structure 

indices 

equipment  support  stiffness 
model  elastic  stiffness 
structure  stiffness  matrix 
equipment  mass 
model  larger  mass 
structure  mass  matrix 
modal  mass 

total  number  of  structural  degrees  of  freedom 
transform  transfer  function  numerator 
Laplace  transform  parameter 
generalized  modal  coordinate 
influence  coefficient  vector 
t ime 

duration  of  input  ground  motion 
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